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This paper summarises an investigation of the statistical properties of orbits escaping from three 
different two-degree-of-freedom Hamiltonian systems which exhibit global stochasticity. Each time- 
independent H = Ho + eH' , with Ho an integrable Hamiltonian and eH' a nonintegrable correction, 
not necessarily small. Despite possessing very different symmetries, ensembles of orbits in all three 
potentials exhibit similar behaviour. For e below a critical eo, escapes are impossible energetically. 
For somewhat higher values, escape is allowed energetically but still many orbits never escape. The 
escape probability P computed for an arbitrary orbit ensemble decays towards zero exponentially. 
At or near a critical value ei > eo there is a rather abrupt qualitative change in behaviour. Above 
ei, P typically exhibits (1) an initial rapid evolution towards a nonzero -Po(e), the value of which is 
independent of the detailed choice of initial conditions, followed by (2) a much slower subsequent 
decay towards zero which, in at least one case, is well fit by a power law P{t) oc t"**, with jj, ~ 0.35 — 

0. 40. In all three cases, Po and the time T required to converge towards Po scale as powers of e — ei, 

1. e., Po oc (e — ei)'* and T oc (e — ei)'', and T also scales in the linear size r of the region sampled for 
initial conditions, i.e., T oc r~^ . To within statistical uncertainties, the best fit values of the critical 
exponents q, /3, and 5 appear to be the same for all three potentials, namely a « 0.5, /3 ~ 0.4, and 
5 ~ 0.1, and satisfy a — (3 — S ~ 0. The transitional behaviour observed near ei is attributed to 
the breakdown of some especially significant KAM tori or cantori. The power law behaviour at late 
times is interpreted as refiecting intrinsic diffusion of chaotic orbits through cantori surrounding 
islands of regular orbits. 



PACS number(s): 03.20.-fi, 05.45.-Fb, 46.10.-Fz 
I. INTRODUCTION 

This paper summarises a study of the problem of es- 
capes of energetically unbound orbits in strongly non- 
integrable two-degree-of-freedom Hamiltonian systems as 
an example of phase space transport in complex sys- 
tems. This work has led to two significant conclusions: 
(1) When evolved into the future, ensembles of orbits 
of fixed energy often exhibit a rapid approach towards 
a constant escape probability Pq, the value of which is 
independent of the details of the ensemble and exhibits 
interesting scaling behaviour. Moreover, the values of 
the critical exponents appear to be relatively insensitive 
to the choice of Hamiltonian. (2) At later times, the es- 
cape probability decreases in a fashion which, for at least 
one model system, is well fit by a power law P{t) cx t^^ 
with /i « 0.35 — 0.40. This nontrivial time-dependence 
is attributed to the fact that the possibility of escape to 
infinity is controlled by cantori, which can trap chaotic 
orbits near regular regions for extremely long times. 



The first three papers in this series (hereafter Pa- 
pers 1-3) described a numerical investigation of the sta- 
tistical properties of orbit ensembles evolving in noninte- 
grable two-degree-of-freedom Hamiltonian systems where 
it is possible energetically for trajectories to escape to in- 
finity. Earlier investigations of individual orbits in these 
systems had led to two significant conclusions [^1): (1) 
Just because escape is possible energetically does not 
mean that escape will inevitably occur; and, even if es- 
cape does occur, the time required for a trajectory to 
cross a Lyapunov curve Q and hence escape to infinity 
may be very long compared with the natural crossing 
time. (2) In some phase space regions, the time and di- 
rection of escape vary smoothly as a function of initial 
conditions, but in other cases one discovers instead an 
apparent near-fractal dependence on the specific choice 
of initial data (cf. Q). 

This problem of escapes is an example of phase space 
transport in complex two-degree-of-freedom Hamiltonian 
systems, a subject which has been explored in detail 
over the past two decades (cf. P-pdj and references cited 
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therein) . The passage of orbits through Lyapunov curves 
and their subsequent escape to infinity is the most con- 
spicuous aspect of the transport, but crucial features of 
the bulk flow, especially at late times, appear to be con- 
trolled by diffusion through cantori |l^-[l5[|, which can 
trap orbits for very long times (cf. |^). The chaotic be- 
haviour of late escapers indicates that this pi'oblem is 
closely related to chaotic scattering (cf. |l6| , |l7|| ), where 
an incident trajectory scatters to infinity at a time and 
in a direction that can exhibit a fractal dependence on the 
value of the impact parameter. However, the problem of 
escapes is also related to a variety of other physical prob- 
lems, including, e.g., the dissociation of molecules (cf. 
10 ) or the evaporation of stars from a cluster (cf. ) . 
Indeed, the phase space interpretation of the escape prob- 
lem suggested in Section IV is completely consistent with 
the detailed phase space description deduced recently for 
escapes in the planar isoceles three-body problem. 

Paper I showed that, for at least one particular Hamil- 
tonian system (the i?i of eq. [1]), the microscopic chaos 
exhibited by individual orbits leads macroscopically to 
bulk regularities: (1) For sufficiently large deviations 
from integrability, localised ensembles of initial condi- 
tions evolve so as to exhibit a rapid approach towards 
a near-constant escape probability Pq , which is indepen- 
dent of the specific choice of initial conditions and which, 
if at all, only changes on a significantly longer time scale. 
(2) The value of Pq scales in terms of an "order parame- 
ter" e — ei. (3) For ensembles that probe a phase space 
region of specified size r, the time T required to converge 
towards Pq also scales in e — ei. (4) For fixed values of 
e — ei , T also scales in the size r of the phase space region 
sampled by the initial ensemble. 

Subsequent work, summarised and extended here, has 
addressed several questions not considered in Paper 1: 

(1) Do other potentials exhibit similar behaviour and, if 
so, are the scaling exponents the same? In other words, 
could this behaviour be universal? 

(2) What happens at much later times? Does P remain 
constant or is there a different asymptotic behaviour for 
larger values of tl 

(3) Is there any obvious correlation between the time 
of escape for individual orbits in the ensemble and the 
exponential instability of those orbits, as probed, e.g., by 
short time Lyapunov exponents (cf. plj-p^)? 

(4) Can one identify a simple, physically well motivated 
model to explain the observed behaviour? 

This work is based on a detailed investigation of orbits 
in three different Hamiltonians, namely: 

Hi = ^{i^+f + x^ + j/)-ex^y\ (1) 
H2^^{x^+f+x^ + y^)-exy^ (2) 



and 

^3 - ^ (^x' + f +x^+y^- ^2/^^ -f ex^y, (3) 

where x and y denote canonical momenta. In each case, 
the Hamiltonian is of the form 

H = Ho + eH', (4) 

with Hq integrable and eH' a nonintegrable correction. 
e is taken to be non-negative but is not assumed to be 
small. These Hamiltonians exhibit very different symme- 
tries. Hi is invariant under x — > —x and/or y — s- —y and 
has four identical escape channels. H2 is only symmetric 
with respect to y — *■ —y, and has two channels of escape. 
For e = 1, H^ reduces to motion in the Henon-Heiles po- 
tential, which manifests a 27r/3 rotation symmetry, but 
for e > 1 this discrete symmetry is broken. Representa- 
tive equipotential surfaces are exhibited in Fig. 1. 

This research has led to four principal conclusions: 
(1) For all three systems, there exists a critical ei, larger 
than eo, the smallest e for which escapes can occur, which 
signals a qualitative change in short time behaviour: Be- 
low ei the escape probability P decays towards zero ex- 
ponentially, but for larger values of e one observes instead 
an initial approach towards a near-constant nonzero Pq, 
the value of which is independent of the detailed choice 
of initial conditions. 

2) For all three systems, Pq scales in e— ei, i.e., Pq oc (e — 
ei)" with a > 0. For a uniform sampling of a given 
phase space region of fixed size r, the time T required to 
converge towards Pq also scales, i.e., T cx (e — ei)^'^, with 
/3 > 0. For fixed e, T also scales in the linear size r of the 
phase space region initially sampled, i.e., T cx r^* with 
S > 0. Finally, the data are at least consistent with the 
possibility that the numerical values of a, /3, and i5 are 
the same for all three potentials; and that a — P — 6 = 0. 

(3) At later times, P deviates from Pq by exhibiting a 
slow decrease towards zero. For at least one Hamiltonian, 
namely H3, this later time evolution is well fit by a power 
law P oc t~^, with a positive constant /i < 1. 

(4) At least for H3, and possibly for Hi and H2, orbits 
that escape early on tend to be more unstable than orbits 
which only escape at much later times, in that they have 
a larger short time Lyapunov exponent. Computed dis- 
tributions of short time Lyapunov exponents and surfaces 
of section both support the hypothesis that the chaotic 
orbits divide, at least approximately, into relatively dis- 
tinct subclasses, presumably separated by cantori. 

The last two points suggest strongly that the late time 
evolution of P(t) is controlled by cantori, which can trap 
chaotic orbits near regular islands for very long times. It 
is well known (cf. [^) that cantori typically constitute 
the principal impediment for efficient phase space trans- 
port in two-degree-of-frecdom systems and that diffusion 
through cantori is usually not characterised by a constant 
escape probability. 
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Section II describes the observed short time behaviour, 
summarising results presented in Papers 1-3 and dis- 
cussing the evidence for universahty. Section III focuses 
on longer time evolution, using surfaces of section and 
short time Lyapunov exponents to provide insights into 
flows associated with initially localised orbit ensembles. 
Section IV suggests a tentative physical interpretation 
of the numerical results in terms of flows in a chaotic 
phase space partitioned by cantori. Phase space is pre- 
sumed to be dominated by "unconfined" chaotic orbits 
which, in the absence of any cantori or Lyapunov curves, 
would evolve towards a statistical equilibrium (cf. 
The constant Pq observed at relatively early times is 
attributed to the fact that orbits sampling this near- 
equilibrium will escape at a near-constant rate. The de- 
caying P{t) later on reflects the fact that the phase space 
also includes an appreciable measure of temporarily "con- 
fined" or "sticky" (cf. | |l5| , p6t ) chaotic orbits which, albeit 
not trapped within the system forever, can only escape 
much later once they have breached one or more cantori 
to become unconfined. 



II. SHORT TIME BEHAVIOUR 
A. Description of the experiments 

The experiments described here entailed a study of or- 
bit ensembles with fixed energy h evolving in Hi, H2, 
and H3 with variable e. Attention focused exclusively on 
e > eo{h), the smallest value of e for which escape to in- 
finity is possible energetically. The values of h and ep for 
all three cases are given in Table 1. 

Ensembles of initial conditions were generated by uni- 
formly sampling a square cell of linear dimension r in 
the {x, x) plane, setting y — 0, and then computing 
y = y{x,i, h) > 0. The cells were chosen to sample "in- 
teresting" phase space regions where most of the orbits 
do not escape at very early times. This implied that the 
time of escape for any given orbit typically manifested a 
sensitive dependence on initial conditions. Most of the 
computations involved a fiducial cell size r — 0.05. How- 
ever, when exploring the eff'ects of varying cell size, cells 
as small as r = 1.0 x 10~^ were also used. 

Each initial condition was integrated into the future 
and the location of the orbit on a Poincare section noted 
at successive consequents, i.e., successive crossings of the 
X = phase space hyperplane with x < 0. If after con- 
sequent t — 1 but before consequent t the orbit crossed 
one of the Lyapunov curves, the orbit was recorded as 
having escaped at t. The experiments with Hi [|l|,§j5| 
used a time series integrator. However, the experiments 
involving H2 and H3 g,^ exploited a more efficient Lie 
integrator truncated at twelfth order |2^], which facili- 
tated integrations of significantly larger ensembles typi- 
cally containing 2000 x 2000 initial conditions or more. 



The fundamental object of interest is P{e,t), the prob- 
ability that a randomly chosen orbit is an escaping orbit 
with escape occuring between consequents t — 1 and t. 
This escape probability, along with an estimated uncer- 
tainty AP{e, t), was defined by the obvious relation 



P(e,i)±AP(e,i) = ^:^^^±^^, (5) 

where iVesc denotes the number of trajectories that es- 
cape between i — 1 and t and Ntot the total number 
present at i — 1 . 

B. Results from the experiments 

At early times, the escape probability can exhibit a 
complex, highly irregular behaviour, the details of which 
depend sensitively on the size and location of the initial 
cell. However, at somewhat later times, P tends instead 
to exhibit a more systematic behaviour that is seemingly 
independent of the initial cell and depends only on the 
value of e. The qualitative form of this behaviour depends 
crucially on whether e is above or below a critical ei > eg. 

Below ei, the escape probability P(e, t) decays towards 
zero in a fashion that is well fit by an exponential. This 
is, e.g., illustrated in Figs. 2 a and b, which exhibit P{t) 
and InP(i) for one value of e < ei in H2, namely e = 
1.04. The solid curve superimposes a best fit exponential 
P cx exp(-i/r), with r = 62.0. 

For e > ei, P{e,t) appears instead to converge to- 
wards a nonzero Po(e)i the value of which depends on 
e but is independent of the cell of initial conditions. The 
evidence is particularly compelling for Hi and H^, but 
somewhat weaker for H2 where, especially for small e— ci, 
the convergence is relatively slow and can merge into the 
later time evolution described in Section III. Examples 
of this behaviour are provided in Figs. 3 a-c, which ex- 
hibit P{t) for representative values e > ei in Hi, H2, and 
H3. (Other examples are provided in Papers 1 - 3.) The 
transition at or near ei is quite abrupt and, for e > ei, 
Po{e) is a monotonically increasing function of e. More- 
over, for all three Hamiltonians one finds that, at least 
for relatively small values of the order parameter e — ei , 
the escape probability Po{e) exhibits a simple scaling, 
namely 

Po(e)«(e-eir, (6) 

with a constant a > 0. Illustrations of the goodness of 
fit of this scaling relation for the Hamiltonians Hi and 
H3 are provided, respectively, by Figs. 7 in Paper 1 and 
Figs. 4 in Paper 3, which exhibit plots of P vs. e and 
InP vs. In (e — ei). The values of ei are again given in Ta- 
ble 1. That the escape probability approaches a roughly 
time-independent value can be interpreted by supposing 
that the cell of initial conditions has dispersed to fill cer- 
tain large regions inside the Lyapunov curves in a nearly 
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uniform fashion, and that trajectories are escaping near- 
randomly at a constant rate. 

Using operational prescriptions described in Paper 1, 
one can also estimate the time T required for an ensemble 
with e > ei to approach Pq. For cells of fixed size r, this T 
appears to be roughly independent of initial conditions, 
depending only on e. Moreover, one finds that, for all 
three Hamiltonians, T also scales in e i.e.. 



T(e)(x(e-ei) 



-0 



(7) 



with /3 > 0. Illustrations of the goodness of fit of this 
scaling relation for the Hamiltonians Hi and are pro- 
vided, respectively, by Figs. 8 in Paper 1 and Figs. 5 
in Paper 3, which exhibit plots of T vs. e and InT 
vs. In (e — ei). 

To the extent that a constant Pq reflects a popula- 
tion that has dispersed throughout the regions inside the 
Lyapunov curves, one might anticipate that the conver- 
gence time T would depend on the size r of the initial cell, 
smaller cells approaching Po(e) more slowly. This too was 
confirmed numerically. Indeed, for all three Hamiltonians 
one finds that, for a fixed value of e, the convergence time 
T(r) also scales in r 



T{r) cx r 



(8) 



with (5 > 0. This is illustrated in Figs. 9 and 10 of Paper 
1 for two different values of e for the Hamiltonian Hi . 

It is obviously important to determine how abrupt the 
transition at e « ei actually is. However, this is diffi- 
cult numerically. Eq. (7) implies that T{e) diverges for 
e — > ei , but this critical slowing down implies an intrinsic 
limitation in one's ability to probe the details near the 
transition point. 



C. Possible evidence for universality 

That the scaling relations (6) - (8) hold for all three 
Hamiltonians is clearly interesting. Even more striking, 
however, is the fact that, as is evident from Table 1, the 
values of the exponents a, (3, and 5 are very similar for 
all three systems. In each case. 



0.5, /3~0.4, 



and 



0.1. 



(9) 



The uncertainties in a are dominated by uncertainties 
in the correct value of ei. As discussed in Paper 1, be- 
cause of the aforementioned critical slowing down ei is 
best estimated by looking at somewhat higher values of 
e and extrapolating to smaller e. The uncertainties in (3 
are dominated by the precise prescription used to identify 
a convergence time. In particular, even though it is usu- 
ally easy to determine a lower bound on the convergence 
time, the determination of an upper bound proves more 
difficult. It follows that the quoted error bars in Table 1 
can be asymmetric. As regards the best fit (5, there are 
two principal sources of uncertainty, namely (1) that the 



effect is relatively small (so that the fractional error is 
large) and (2) that, especially for very small values of r, 
different ensembles can exhibit significant variability. 

Given these uncertainties, one cannot conclude unam- 
biguously that a, f3, and S are strictly equal for all three 
Hamiltonians. However, one can conclude that they are 
all comparable in size for all three systems and that, con- 
sistent with the uncertainties, they may in fact be equal. 

It is also true that, to within statistical uncertainties. 



a- (3- 6^0. 



(10) 



For Hi and H^ the evidence for this assertion is relatively 
strong. For H2 the case is somewhat weaker, largely be- 
cause the best fit /3 is somewhat larger than for the other 
two systems. However, it should be noted that the error 
bars on the f] for H2 are especially big. This reflects the 
fact that, for this system, the short time behaviour de- 
scribed in this Section merges relatively quickly into the 
later time evolution described in Section HI, where P{t) 
begins to decay to values below Pq. For Hi and H3 this 
subsequent decay only become significant at somewhat 
later times, at least for larger values of e — ei. 

The evidence for universality summarised here is cer- 
tainly much weaker than for the universality first identi- 
fied by Feigenbaum |28| for one-dimensional maps or by 
Escande and Doveil 1 29 1 and MacKay |Q in their renor- 
malisation group analyses of tori, but it is, nevertheless, 
intriguing. Moreover, certain points are seemingly un- 
ambiguous. (1) For all three Hamiltonians, there is clear 
evidence for an abrupt change in behaviour at or near 
some critical value ei. (2) For e > ei, orbit ensembles 
evolve towards an escape probability P{t) which is (a) 
largely independent of the choice of initial ensemble and 
(b) nearly time- independent, at least for relatively short 
times. (3) The convergence time T depends both on cell 
size r and e — ei, larger cells and larger e — ei leading to 
a more rapid convergence. (4) Because T(e, r) increases 
with decreasing e — ei, pinning down the precise value 
of ei is quite hard. (5) The resulting uncertainties in ei 
do not impact the apparent fact that P and T scale in 
e — ei. However, they do impact estimates of the precise 
values of a, f3, and 5, thus making it difficult to deter- 
mine whether these exponents are the same for all three 
potentials. What is clear is that, for all three potentials, 
the values of the exponents are comparable in magnitude. 



III. PHASE SPACE FLOW AT LATER TIMES 

A. Late time evolution of P{t) 

Section II summarised experiments indicating that, on 
a relatively short time scale, the escape probability P{t) 
evolves towards a near-constant value Poi^)- However, 
there is no reason to expect that P will remain approx- 
imately constant if the orbit ensembles are evolved for 
much longer times. If, e.g., the initial ensembles contain 
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a few regular orbits that cannot escape, these will even- 
tually dominate the orbits that remain in the system and 
P{t) must decay towards zero exponentially. Even if the 
initial ensembles contain no regular orbits that never es- 
cape, one might anticipate more complicated behaviour 
at late times. If, e.g., some small subset of the chaotic es- 
cape orbits are stuck by cantori near some regular island 
for relatively long times, the escape probability should 
decrease below the initial near-constant Pq once the other 
chaotic orbits have almost all escaped. 

Such a decrease was first noted for orbit ensembles in 
H2 and subsequently studied more systematically for en- 
sembles in H3 |3|. The principal conclusion is that, at 
least for e > ei, when orbit ensembles are evolved for 
somewhat longer times the probability P{t) begins to ex- 
hibit a slow, monotonic decrease. This is illustrated in 
Fig. 4 which, for one ensemble evolved in H^, exhibits 
P{t) for an ensemble with '--^ 1.1 x 10^ orbits carefully 
chosen from a tiny region of size r — Ix 10^^ to be dom- 
inated by "slow escapers," so that P{t) can be tracked for 
a comparatively long interval. Visually, P{t) decays too 
slowly, and has the wrong curvature, to be well fit by an 
exponential. However, as is illustrated in Fig. 4b, which 
plots In P as a function of \nt, the data for 5 < t < 180 
or so can be fit quite well by a power law. 



P{t) cx r 



(11) 



with fi = 0.39±0.02. 

This algebraic decay seems very robust, with fi appar- 
ently independent of the initial cell and, at least within 
a limited range, the value of e. The best fit value for 
several different values of e (most of which were sampled 
for shorter times with far fewer orbits) is /i = 0.35±0.07. 



B. Tools of analysis 

To ascertain why P{t) changes in time and to better 
understand the qualitative character of the flow, orbit 
ensembles evolved in H3 were also analysed in two other 
ways. 

The first involved computing surfaces of section for an 
evolving ensemble. Each orbit in the ensemble was inte- 
grated into the future and, provided that it had not yet 
escaped, its values of y and y were recorded at successive 
consequents and sorted to generate sequences of surfaces 
of section exhibiting (y, y) pairs. 

Such surfaces of section allow one to determine the ex- 
tent to which the ensemble has evolved to cover a large 
portion of the allowed phase space. Moreover, they can 
facilitate the detection of "zones of avoidance" associated 
with regular islands and/or with orbits that have imme- 
diately escaped, as well as phase space regions with ex- 
cess concentrations of orbits, corresponding presumably 
to regions from which escape is especially difficult. 

The second involved computing short time Lyapunov 
exponents (cf. pl|-p3|), which probe the average expo- 



nential instability of chaotic orbits over finite time inter- 
vals. By analogy with ordinary Lyapunov exponents, a 
finite time can be defined by the obvious prescription 



yU) = lim — In j ^ ^ f \ I 
A-Z(0)^0 t \SZ{0)\ 



(12) 



where = {Sx)^ + {5yf + {5xY + {5yY denotes 

the magnitude of the initial phase space perturbation, 
defined with respect to the natural Euclidean norm. 
These exponents were determined computationally in the 
usual way [p5| by introducing a small initial perturbation 
(5x(0) = 1 X 10~^^ in the x-direction and evolving simulta- 
neously both the perturbed and unperturbed orbits, peri- 
odically renormalising the amplitude of the perturbation 
to assure that |52'(t)| remains smaller than 1 x 10^^. 

At early times, the x{t) computed in this way will de- 
pend strongly on the initial perturbation. However, if 
the orbit segment is chaotic, x(t) will quickly become 
dominated by the component of the perturbation in the 
most unstable direction and become relatively insensitive 
to the initial 5Z{Q). Note also that, given x(i) for two 
different times, ti and ^2, one can identify the average 
exponential instability for the interval ti < t < t2 ss 



X{t2-tl) 



hxih) - tixjti) 



(13) 



Short time Lyapunov exponents were used to confirm 
that, for the values of h and e under consideration, most, 
if not all, of the computed orbits in are chaotic. Dis- 
tributions of short time Lyapunov exponents were also 
used to show (1) that the chaotic orbits which have not 
escaped often appear to divide into distinct populations 
and (2) that there are correlations between the magni- 
tude of X and the time at which the orbit escapes. 



C. Surfaces of section 

Figures 5 a-f exhibit a sequence of sections, generated 
for one ensemble at six different consequents, t = 2, 6, 10, 
15, 20, and 25. This ensemble was comprised of 160, 000 
orbits evolved with e = 1.13, a value only marginally 
above the critical ei. The cell of initial conditions, with 
< x < 0.05, 0.04 <x< 0.09, y = 0,h= 1/6, and hence 
0.5681 < y < 0.5780 is located near the top of the ener- 
getically accessible portion of the surface of section. The 
first escapes occured at t = 4, and the largest number of 
escapes was at i = 7. P{t) first settled down towards a 
smoothly varying form around f = 10 — 12. 

Inspection of these, and other intermediate, sections 
indicates that the orbits remaining inside the Lyapunov 
curves tend systematically to spread over a relatively 
large fraction of the energetically allowed phase space. 
Indeed, the elongated striae associated with the specific 
choice of initial conditions, so conspicuous at consequents 
t — 2 and 6, have been significantly blurred by t = 10 
and have nearly disappeared by i = 25. 
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For t < 6 — 8, the form of the sections is strongly time- 
dependent. However, by t = 10 the ensemble has evolved 
to yield sections characterised by three seemingly dis- 
tinct regions which persist to later times, namely: (1) two 
large holes accompanied by smaller surrounding whorls, 

(2) several overdense regions at positive values of y, and 

(3) a larger region characterised by a substantially lower 
density. As time passes, the occupied regions all decrease 
in density, but the overdense regions remain overdense. 

The two holes and their surrounding whorls are asso- 
ciated with escapes through Lyapunov curves: any orbit 
with values of y and y in these regions would already 
have escaped before intersecting the x-axis. For larger 
values of e it is apparent that the visible whorls are part 
of an elaborate set of structures that penetrate through- 
out large portions of the lower density regions, and that 
this lower density region, which appears macroscopically 
to be populated in a near-uniform fashion, is really laced 
with tiny zones of avoidance. 

Assuming that essentially all the orbits in the ensem- 
ble are chaotic, the overdense regions can be interpreted 
as reflecting trapping near regular islands: Even though 
these islands may be so small as to be almost unobserv- 
able, cantori can significantly impact relatively large por- 
tions of the chaotic phase space HJlltl. The idea then is 
that orbits in the overdense regions are trapped by can- 
tori and can only escape once they have diffused through 
the cantori and can travel unimpeded throughout the re- 
mainder of the chaotic sea. 

The existence of this three-part structure ~ holes, less 
dense regions, and more dense regions - is independent 
of the choice of initial conditions. Moreover, the gen- 
eral locations of the holes and the higher density regions 
are insensitive to the precise value of e. This latter fact 
is manifested in Fig. 6, which shows the analogue of 
Fig. 5 c, now generated for an ensemble with e = 1.06. 
This reinforces the interpretation that one is seeing the 
effects of basic phase space structures, rather than tran- 
sient streaming motions reflecting the choice of ensemble. 



D. Short time Lyapunov exponents 

Perhaps the most obvious way to search for correla- 
tions between the degree of exponential instability exhib- 
ited by a chaotic orbit and the time at which it escapes 
is to compute the mean short time Lyapunov exponent, 
(x(^-e))j for all the orbits in an ensemble that escape be- 
tween successive consequents ts and ts + The results 
of one such computation are presented in Fig. 7, which 
was generated from an ensemble of 1 x 10^ orbits evolved 
in H3 with e = 1.30. {x{tE)) clearly exhibits an initial 
relatively rapid decrease for t^; < 5 — 10, followed by 
a more extended period during which (x(i_B)) decreases 
more slowly. (The initial point in Fig. 7 at = 5 is 
statistically significant.) 

As described above, interpreting the computed (xC^-e)) 



at very early times as an accurate probe of the average 
maximum short time exponent is suspect. Given, how- 
ever, that the typical values of x are greater than or of 
order unity, the computed {xits)) should be relatively re- 
liable for Ie > 5 or so, which means that the initial rapid 
decrease is most likely real. It is not completely clear 
whether will continue to decrease at late times, 

or whether it asymptotes towards a nonzero value. How- 
ever, the existence of a continued decrease out to at least 
t = 70 or so is unquestionably significant statistically. 

The observed decrease in {x{tE)) can be easily inter- 
preted by assuming that the chaotic phase space inside 
the Lyapunov curves divides into two different regions, 
namely (1) a region where orbit segments are less unsta- 
ble exponentially and from which direct escape to infinity 
is difficult, if not impossible, and (2) a region where or- 
bit segments are more unstable and from which escape 
to infinity can proceed on a relatively short time scale. 
The idea is that orbits which remain inside the Lyapunov 
curves for a long time will typically spend most of their 
time in the less unstable region before entering the more 
unstable region and subsequently escaping to infinity. It 
follows that, for orbits that only escape at late times, the 
short time exponent x(^-e), which probes the average in- 
stability for < t < ts, will typically be smaller than for 
early escapers. 

Suppose, oversimplistically, that orbit segments in the 
less and more unstable regions can be characterised re- 
spectively by unique short time exponents xl and xhi 
and that any orbit that enters the high x region will es- 
cape after exactly a time ig- It then follows from eq. 
(13) that the total x(^b) for an orbit escaping at time Ie 
satisfies 



X{tE) = ^ 
tE 



te XH + {tE - te)XL 



(14) 



Consistent with Fig. 7, this implies that x(^-e) de- 
creases monotonically, but eventually asymptotes to- 
wards a nonzero xl- For t^; > 15 or so the computed 
(x(^b)) exhibited in Fig. 7 is in fact well fit by eq. (14) 
with XL = 0.92 and te{xH - Xl) = 7.80. 

To confirm that the observed decrease in (x(^£;)) re- 
flects transitions between relatively distinct orbit popu- 
lations, it is also instructive to determine how short time 
Lyapunov exponents computed for the same set of orbits 
change as a function of time. One way to do this is to 
consider all the orbits in an initial ensemble that escape 
at a given time Te and, given expressions for x(i) for 
times t < Te, computed using eq. (12), extract short 
time exponents for different intervals ti < t <t2. 

Figure 8, generated from the same ensemble as Fig. 7, 
focuses on all 832 orbits that escaped at t = 50, comput- 
ing distributions of short time exponents for the intervals 
15 < i < 20 (solid curve) and 45 < t < 50 (dashed curve). 
Both distributions are bimodal, seemingly comprised of 
two different populations with peaks at or near the same 
values of x, but it is clear that the relative height of the 
two peaks changes significantly in time. For the earlier 
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interval, the lower x population dominates whereas for 
the later interval the low and high x populations are of 
more nearly equal importance. For 10 < i < 15 the mean 
(x) = 0.70; for 45 < t < 50 the mean (x) = 1.11. This is 
consistent with the interpretation that most of the orbits 
that escaped at i = 50 were members of a low x popula- 
tion at early times but shifted to the higher x population 
shortly before escaping from the system. 

IV. PHYSICAL INTERPRETATION 
A. General considerations 

This section suggests a simple model for the behaviour 
observed in Sections II and III which is based on the 
assumption that, over finite intervals, orbits divide at 
least approximately into three distinct classes, namely 
(1) regular orbits, (2) sticky, or temporarily confined, 
orbits, and (3) unconfined chaotic orbits, even though 
the distinction between confined and unconfined chaos 
disappears entirely in the t oo limit. 

Conventional wisdom would suggest that the presence 
of stable periodic orbits, which one expects for generic 
Hamiltonians, implies that there must exist a finite mea- 
sure of regular orbits, even at very high energies h and/or 
large values of e. However, the existence of such regular 
regions, separated from the remaining chaotic orbits by 
invariant KAM tori, suggests in turn that the surround- 
ing chaotic sea should contain cantori which, albeit not 
impenetrable barriers, can trap chaotic orbits near the 
regular islands for relatively long intervals of time. Even 
if there is no absolute distinction between different types 
of orbits in the chaotic sea, there may exist short time de 
facto distinctions which can have significant implications 
for the Hamiltonian flow. (Strictly speaking, in general 
there will also exist a finite measure of chaotic orbits in- 
side the KAM toii. However, these can never breach the 
invariant tori and, as such, will never be able to escape 
to enter the surrounding stochastic sea. For this reason, 
they may be lumped together with the regular orbits in 
the following discussion.) 

If a localised ensemble of phase space points, each cor- 
responding to a chaotic orbit with energy h, is evolved 
in a time-independent potential which, as for the Hamil- 
tonians (1) - (3) for e < eoi has a compact constant en- 
ergy hypersurface, one anticipates a coarse-grained evo- 
lution towards an invariant distribution corresponding 
to a uniform population of the accessible phase space 
[ pT| . Numerical experiments suggest (cf. [^) that, if 
the chaotic phase space is not significantly impacted by 
cantori, so that a single chaotic orbit can easily access 
the entire region without having to diffuse through any 
barriers (cf. p2| , ^ ), this approach will proceed expo- 
nentially in time on a time scale of order the natural 
crossing time, tcr- If, however, cantori play an important 
role in partitioning the chaotic phase space regions over 



relatively short time scales, one can instead observe a 
more complex, seemingly two stage process [^-p4| . Sets 
of orbits in different nearly disjoint regions will rapidly 
approach near-invariant distributions, corresponding to 
near-uniform populations of the separate regions; but 
only later, on a significantly longer time scale, will or- 
bits from different regions "mix" to yield an approach 
towards a true invariant distribution. 

Suppose now that e > eo, so that orbits are no longer 
bound energetically. It then seems reasonable to inter- 
pret the observed behaviour of orbits in terms of two 
distinct sorts of "escape," namely (1) unconfined chaotic 
orbits which pass through Lyapunov curves to escape to 
infinity and (2) confined chaotic orbits, originally stuck 
near the regular regions, which pass through one or more 
cantori to become unconfined, after which they too can 
escape to infinity. To the extent that the escape chan- 
nels - both the simple gaps breached by Lyapunov curves 
and the more complex cantor sets of holes in cantori - are 
small, or that one is considering phase space regions rel- 
atively far from the escape channels, it should be reason- 
able to visualise what is happening in terms of ensembles 
that have evolved towards a near-invariant distribution. 

Suppose, in particular, that one selects an initial en- 
semble where most of the orbits are unconfined chaotic, 
but that there are also a significant number of confined 
chaotic orbits and, perhaps, a few regular orbits. It is 
then easy to explain the qualitative evolution. 

On a relatively short time scale, the unconfined chaotic 
orbits should evolve towards a near-uniform population 
of the phase space regions far from the Lyapunov curves 
and outside any cantori which significantly impede phase 
space transport. (The fact that, in the late time sections 
of Fig. 5, different parts of the lower density region seem 
to have the same relative density at different times sup- 
ports this idea.) However, once this near- invariant dis- 
tribution has been achieved, escapes to infinity should 
proceed "at random" at a near-constant rate, so that the 
unconfined chaotic orbits will be characterised by a con- 
stant escape probability. To the extent that the total 
orbit population is dominated by the initially unconfined 
orbits, and that appreciable numbers of confined orbits 
have not yet become unconfined, the total escape proba- 
bility P should be approximately constant. 

Eventually, however, most of the unconfined chaotic 
orbits will have escaped, so that, assuming that uncon- 
fined orbits cross the Lyapunov curves more quickly than 
confined orbits diffuse through cantori, the total escape 
probability is impacted significantly, and ultimately dom- 
inated, by the remaining orbits. If all these orbits were 
regular and unable to escape, one would expect P to de- 
cay exponentially in time (cf. |^]). Given, however, that 
most of the remaining orbits are chaotic orbits originally 
trapped by cantori, as seems true for e > ei , one expects 
a slower decay in P reflecting transitions from confined 
to unconfined chaos: The idea here is that confined or- 
bits will only diffuse through cantori to become uncon- 
fined very slowly, but that, once unconfined, they will 
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quickly escape to infinity. The observed P cx t~'^ is thus 
driven by the diffusion of orbits through cantori rather 
than their escape through Lyapunov curves. That diffu- 
sion through cantori need not characterised by a constant 
escape probabihty is in fact weU known (cf. |9j). 

The same quahtative picture should remain valid even 
if the initial orbit ensemble is carefully selected to contain 
no chaotic orbits trapped near the regular regions. Even 
if most of the original unconfined orbits quickly escape 
to infinity through one of the Lyapunov curves, a small 
fraction of those orbits could become trapped by cantori. 
However, once trapped most of these orbits will only leak 
out at significantly later times, when the population of 
unconfined orbits inside the Lyapunov curves has become 
significantly reduced. 



B. A simple model 



Similarly, Nu{t) must satisfy a simple rate equation 



dt 



dt 



dN^{t) 



dt 



-XN^{t)~NM 



dm 

dt 



This latter equation is easily solved to yield 



Rit) 



(Nr + Nc + Nu) dt 



(17) 



Nu{t) = Nu{0)e-^' - Nc{0) f dre^f^-*) ^ . (18) 

However, this implies that the total escape rate 
-1 d 



Consider an ensemble of initial conditions of energy 
h, located inside the Lyapunov curves, which may be 
divided into three different classes - unconfined chaotic, 
temporarily confined chaotic, and regular - characterised 
by numbers iV„(0) > iVc(O) > Nr{0). Now implement a 
probabilistic description, treating the regular orbits as a 
separate population that can never escape, but allowing 
for three sorts of transitions, namely escapes of uncon- 
fined orbits through Lyapunov curves, trapping of uncon- 
fined orbits by cantori, and leakage of confined orbits to 
become unconfined orbits inside the Lyapunov curves. 

Calculating the total escape rate R, the continuum 
limit of the discrete probability P computed in the nu- 
merical experiments summarised above, is straightfor- 
ward if one makes two basic assumptions, each implicit 
in the preceding and seemingly consistent with the nu- 
merical experiments. (1) The rate at which unconfined 
orbits escape to infinity assumes a constant value A, in- 
dependent of time. (2) A is much larger than the rates 
at which unconfined orbits become temporarily confined 
and temporarily confined orbits become unconfined. 

These assumptions imply that, at early times, the total 
escape rate is dominated by the escape of initially uncon- 
fined orbits, and that the details of any early trapping of 
unconfined orbits are irrelevant. This means that, when 
computing the total escape rate, the confined population 
may be approximated by an expression of the form 



N,{t) = N,{0)f{t), 



(15) 



where Nc{0) allows for a possible early time trapping of 
some small fraction of the originally unconfined orbits 
and f{t) is a monotonically decreasing function of time, 
assumed to satisfy an initial condition /(O) = 1. The 
rate Rcu at which temporarily confined orbits become 
unconfined is thus 



Rcu{t) 



Idl 
fdt 



(16). 



A 


1 - V, /o dTe^^df{T)/dT 




1 - i^c Jo dTe^^df{T)/dT + Vcf{t)e^^ + z^^e^* 



(19) 

where Vc = Nc{G)/Nuifi) and Vr = Nr/Nu{Q) reflect orig- 
inal relative abundances. 

That the escape rate through the Lyapunov curves 
in much larger than the rate at which confined or- 
bits become unconfined means that A ^ \(\/ f){df / dt)\. 
This implies, however, that the integrals in the preced- 
ing equation can be evaluated perturbatively, expanding 
df{T)/dT about its value at time t = t. Recognising that 
the second term in the denominator is small compared 
with the third, one thus concludes that 



R{t) 



A- 



df{t)/dt 



1 + lycc^'fit) + 



(20) 



Granted that I'r ^ i^c ^ 1, R{t) typically admits three 
different, asymptotic regimes, namely 

(1) an early time regime where R ~ X; 

(2) an intermediate regime where R —{l/f)df/dt; and 

(3) a late time regime where R —{vc/i^r)df/dt. 

At early times, the total escape rate is fixed by the rate 
at which initially unconfined orbits cross the Lyapunov 
curves. Later, once most of these original unconfined 
orbits have escaped, R{t) is set by the rate at which con- 
fined chaotic orbits become unconfined. Finally, once 
most of the initially confined orbits have escaped, there 
is a more rapid decrease in R towards zero, reflecting the 
fact that the now dominant regular population can never 
escape. 

It is natural to interpret the probability P cx 
described in Section HI as reflecting the intermediate 
regime. If /i were equal to unity, one would then infer 
a population 



Nc{t) (X r 



(21) 
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for some constant q > 0, i.e., an eventual power law decay 
in the number of confined orbits. Given, however, that 
the best fit /i < 1, one infers instead a population 



Nc{t) cx exp{~qt 



(22) 



which decays faster than an a power law but slower than 
the exponential decrease associated with a constant es- 
cape rate. 

That P(t) is not constant is hardly surprising. Indeed, 
diffusion through cantori, interpreted as orbits wend- 
ing their way through a self-similar collection of turn- 
stiles, p2|, 33| would suggest that P{t) decay in time 



(cf. I^-^I). Rather, what is interesting is that n 1- 
In the context of chaotic scattering, one seems to see 
a sharp distinction between hyperbolic scattering (cf. 
1^,^), where the number of incident particles remain- 
ing within the scattering region at time t decays expo- 
nentially, i.e., N{t) cx exp(— A<), and nonhyperbolic scat- 
tering (cf. ||39|), where the number remaining decays as 
a power law, i.e., N{t) cx t~'^. The origin of the interme- 
diate behaviour observed here is, at the present, unclear. 

This picture presupposes that ^ J^c ^ 1. Indeed, 
if these inequalities fail, the qualitative evolution can 
change significantly. Suppose, for instance, that tem- 
porarily confined chaotic orbits are unimportant com- 
pared with regular orbits, so that one need consider only 
two orbit classes - regular and unconfined ~, and that, 
even early on, there are many more regular than uncon- 
fined orbits. It then follows immediately that, already 
at early times, R{t) should decay towards zero exponen- 
tially: 



N, 



1 + {Nr/Nc)exp{-Xt) Nr 



exp(-At). (23) 



This is consistent with the observed behaviour for eo < 
e < £i where, as noted already, P{t) decays towards 
zero without first asymptoting towards a near constant 
nonzero value. 



C. Discussion 

To place the preceding in an appropriate context, three 
important points should be noted. 

1. Extracting a near-constant Po{e) is necessarily a some- 
what imprecise operation. The algorithm described in 
Paper 1, or any obvious alternative, depends crucially 
on the idea that P{t) will evolve towards a form that 
is largely independent of the initial conditions on a time 
scale sufficiently short that appreciable numbers of sticky 
chaotic orbits do not escape and the total escape prob- 
ability is dominated by the near-constant rate at which 
unconfined orbits escape. If the ensemble "forgets" its 
initial conditions sufficiently quickly, it is possible op- 
erationally to identify a reasonable estimate of Poie). 
Strictly speaking, however, the ensemble is really evolv- 
ing towards a characteristic P{e, t) which manifests a 



slow, but nontrivial, time-dependence. This leads to an 
inherent inaccuracy in the determination of the values of 
the critical ei and, especially, the exponents a, /?, and S. 

2. The calculations described in Section III were re- 
stricted to relatively short times, t < 200 or so; and, 
for this reason, one cannot preclude the possibility that, 
on a significantly longer time scale, P{t) could assume a 
form very different from what was suggested in Section 
IV. For example, one cannot preclude the possibility that, 
for much later times, P(t) cx , in agreement with Kar- 
ney's | |37[ | experiments. Ideally one might like to integrate 
for much longer times, but this quickly becomes very ex- 
pensive computationally: because the overwhelming ma- 
jority of the orbits in an initial ensemble escape relatively 
early on, one would need to start with an absolutely enor- 
mous collection of initial conditions in order to derive 
statistically significant conclusions about behaviour at 
much later times. Nevertheless, even though one cannot 
exclude the possibility of different later time behaviour, 
it is significant that the observed scaling P{t) cx for 
t < 180 or so appears to be robust 

It should also be noted that one cannot completely ex- 
clude the possibility that the long time integrations de- 
scribed in Section III are contaminated by accumulating 
errors which, e.g., make the system slightly dissipative. 
Numerical tests described in Paper I allow one to be con- 
fident that the early time computations (t < 20 or so) are 
reliable, but there is less hard evidence to justify accept- 
ing long time integrations at face value. All that can be 
said definitively is that, even for the longest time inte- 
grations, the relative energy error for an orbit was never 
larger than 2 x 10~^ and usually orders of magnitude 
smaller. 

3. This simple three-component model, based on two, 
and only two, nearly distinct classes of chaotic orbits, 
may well be oversimplistic. Indeed, lumping together 
every chaotic orbit that is not unconfined into a single 
population assumed to have reached a statistical near- 
equilibrium seems less justified than assuming that, at 
least away from the Lyapunov curves, the unconfined 
phase space is characterised by a near-invariant distribu- 
tion. Furthermore, the entire analysis assumes an abrupt 
transition occuring at or near some critical ei. However, 
one can argue that these limitations are not completely 
unreasonable. Detailed examinations of orbit ensembles 
on a compact phase space hypersurface indicate that, of- 
tentimes, many of the qualitative features of a flow can be 
interpreted by allowing only for two classes of chaotic or- 
bits, named sticky and unconfined (cf. ||3^ , |3^ , ^ ) . More- 
over, investigations of the effects of increasing deviations 
from integrability suggest that, as the control parameter 
e becomes larger, holes in cantori can abruptly increase 
in size, so that what was initially a barrier that could 
only be penetrated on a very long time scale ceases to 
jlay a significant role in impeding phase space transport 
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TABLE I. Critical escape parameters 




FIG. 1. Equipotential surfaces of the potential for (a) Hi 
with £ = 5.26, (b) H2 with e = 3.0, and (c) Ha with e = 1.0. 
Oi, O2, O3 and O4 represent Lyapunov unstable periodic or- 
bits. 
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FIG. 2. (a) The escape probability, P{t), computed for 
five different cells of initial conditions evolved in H2 with 
e = 1.04 < ei. The solid curve is an exponential fitted to 
the interval 20 < i < 100. The error bars reflect uncertainties 
computed as in eq. (5). (b) The same data on a log-log plot. 
Here the orbits from the different cells have been combined 
to yield a single P and AP, again computed as in eq. (5). 
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FIG. 3. (a) The escape probability, computed 
for several different cells of initial conditions in Hi with 
e = 5.05 > ei. The dashed line indicates the best fit Pq. 
(b) The same for H2 with e = 1.20. (c) The same for H3 with 
e = 1.30 
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FIG. 4. (a) The escape probability, P{t), computed for 
the same initial conditions as in Fig. 3 (c), now allowing for 
much longer times. The solid curve exhibits a power law fit, 
P oc with /X = 0.39. (b) The same data on a semilog plot, 
analysed as for Fig. 2b. 




FIG. 5. (a) Surface of section at t = 2 for an ensemble 
of initial conditions evolved in H3 with e = 1.13, exhibiting 
y and y for every orbit still inside the Lyapunov curves that 
crosses the a; = axis with a; < between consequents t and 
t+1. 




(b) The same for i = 6. 





(d) The same for t = 15. 




(e) The same for t = 20. 




(c) The same for t = 10. 



(f) The same for t = 25. 
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FIG. 6. A surface of section analogous to Fig 5 c, generated 
at i = 10 for an ensemble of orbits with e = 1.06. 
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FIG. 8. Distributions of short time Lyapunov exponents 
for the intervals 15 < i < 20 (solid curve) and 45 < i < 50 
(dashed curve) for orbits that escape through the Lyapunov 
curves at consequent t = 50. 



FIG. 7. (a) The escape probability, P{t), computed for 
the same initial conditions as in Fig. 3 (c), now allowing for 
much longer times. The solid curve exhibits a power law fit, 
P oc with jj, = 0.37. (b) The same data on a semilog plot, 
analysed as for Fig. 2b. 



